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Abstract 

Influence of asymmetry on superfluidity of nuclear matter with triplet-singlet 
pairing of nucleons (in spin and isospin spaces) is considered within the frame- 
work of a Fermi-liquid theory. Solutions of self-consistent equations for the 
critical temperature and the energy gap at T = are obtained with the use 
of Skyrme effective nucleon interaction. It is shown, that if the chemical po- 
tentials of protons and neutrons are determined in the approximation of ideal 
Fermi-gas, then the energy gap for some values of density and asymmetry 
parameter of nuclear matter demonstrates double-valued behavior. However, 
with account for the feedback of pairing correlations through the normal dis- 
tribution functions of nucleons two-valued behavior of the energy gap turns 
into universal one- valued behavior. At T = superfluidity arises and disap- 
pears as a result of a flrst order phase transition in density. 
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It is well established, that neutron-proton (np) pairing plays an essential role in the de- 
scription of superfluidity of finite nuclei with N = Z (see Refs. |1|,0] and references therein) 
and symmetric nuclear matter 0-1^. In astrophysical context np pairing correlations can 
be important for the description of r-process and cooling of neutron stars, which per- 
mit pion or kaon condensation In this Rapid Communication we shall investigate the 
infiuence of asymmetry on np superfiuidity of nuclear matter. Previously this problem was 
treated with the use of various approaches and potentials of NN interaction. In particular, 
the cases of ^Si-^Di and pairing were considered in Refs. [0,|T3 on the basis of the 
Thouless criterion for the thermodynamic T matrix. As a potential of NN interaction, the 
Graz II and Paris potentials were chosen, respectively. Superfiuidity in ^Si-^Di pairing 
channel was studied also in Ref. [jll| within BCS theory of superconductivity with the use 



of the Paris potential in the separable form. Investigations, based on the Thouless crite- 
rion, deduce the suppression of np pairing correlations with increase of isospin asymmetry. 
However, the Thouless criterion can be exploited for finding the critical temperature only, 
but does not permit one to draw any conclusions about superfiuidity with a finite gap. The 
studies in Ref. [0, based on the BCS theory, were carried out with the use of the bare in- 



teraction and the single particle spectrum of a free Fermi gas and give, thus, overestimated 

on 



values of the energy gap. The effect of ladder-renormalized single particle spectrum |12 
the magnitude of the energy gap in ^Si-^Di pairing channel was investigated in Ref. |T3| . 
The Argonne V14 potential was explored as input for determination of the single particle 
energy and the bare interaction in the form of the Paris potential was used to evaluate the 
energy gap. The use of the bare interaction in the gap equation seems to be a very strong 
simplification, because medium polarization strongly reduces the magnitude of the gap (see 



Ref. |T^ for the infiuence of the polarization effects on the pairing force in the ^5*0 channel). 
In principle, the effective pairing interaction should be obtained by means of Brueckner 
renormalization, which gives the correct interaction after modifying the bare interaction for 
the effect of nuclear medium. However, the issue of a microscopic many-body calculation 
of the effective pairing potential is a complex one and still is not solved. For this reason, 
it is quite a natural step to develop some kind of a phenomenological theory, where instead 
of microscopical calculation of the pairing interaction one exploits the phenomenological ef- 
fective interaction. We shall investigate the influence of asymmetry on superfluid properties 
of nuclear matter, using Landau's semiphenomenological theory of a Fermi-liquid (FL). In 
the Fermi-liquid model the normal and anomalous FL interaction amplitudes are taken into 
account on an equal footing. This will allow us to consider consistently the influence of the 
FL amplitudes on superfluid properties of nuclear matter. Besides, as a potential of NN 
interaction we choose the Skyrme effective forces, describing the interaction of two nucleons 
in the presence of nucleon medium. The Skyrme forces are widely used in the description 
of nuclear system properties and, in particular, they were exploited for the description of 
superfluid properties of flnite nuclei as well as inflnite symmetric nuclear matter ||T7 



The basic formalism is laid out in more detail in Ref. [0, where superfluidity of sym- 
metric nuclear matter was studied. As shown there, superfluidity with triplet-singlet (TS) 
pairing of nucleons (total spin 5* and isospin T of a pair are equal S = 1, T = 0) is real- 
ized near the saturation density in symmetric nuclear matter with the Skyrme interaction. 
Therefore, we shall study further the influence of asymmetry on superfluid properties of 
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TS phase of nuclear matter. For the states with the projections of total spin and isospin 
Sz = Tz = the normal distribution function / and the anomalous distribution function g 
have the form 



/(P) = /oo(p)c^07-0 + /03(p)cr0T-3, gip) = 5-30 (p)c^3Cr2T-2 (l) 

where cTj, are the Pauli matrices in spin and isospin spaces. The operator of quasiparticle 
energy e and the matrix order parameter A of the system for the energy functional, being 
invariant with respect to rotations in spin and isospin spaces, have the analogous structure 

£(P) = £^00(P)0-07-0 + £^03(P)0^01"3, A(p) = A3o(p)(T30-2r2. (2) 

Using the minimum principle of the thermodynamic potential and procedure of block diago- 
nalization one can express evidently the distribution functions /oo, /03 and 5^30 in terms 



of the quantities e and A: 

^ ^ ~ 4I (^^"^ ^ 2"/°^ ^ ^^^^ ^ 2T ) ' 
/o3 = - ^ (tanh - tanh ^^^) , (4) 

930 = ^ftanh — + tanh zr^V (5) 

y3o AE ^ 2T 2T ^ ^ ^ 

Here 

E = \/^00 + ^30 5 Coo = £^00 - Atoo; Co3 = ^03 - /^035 

T is temperature, Hp and /i° are chemical potentials of protons and neutrons. To obtain the 
closed system of equations for the quasiparticle energy e and the energy gap A, it is necessary 
to express the quantities e, A through the distribution functions / and g. For this purpose 
one has to set the energy functional S{f, g) of the system. In the case of asymmetric nuclear 
matter with TS pairing of nucleons the energy functional is characterized by two normal 
Uo, U2 and one anomalous Vi FL amplitudes ||19[. Differentiating the functional S{f, g) with 



respect to g and using Eq. (^, one can obtain the gap equation in the form 



2T 2T 



The anomalous interaction amplitude Vi in the Skyrme model reads [19 

Vi{p,q) =to{l+Xo) + lt3(f{l + X3) (7) 
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where g is density of nuclear matter, ti,Xi,(3 are some phenomenological parameters, which 
differ for various versions of the Skyrme forces (later we shall use the SkP potential [0). 
Equation (||) should be solved jointly with equations 



?oo(p) 



V^t E{p) 



being the normalization conditions for the normal distribution functions /oo, /os- In Eq- (0) 
the quantity a = {gn — gp)/g is the asymmetry parameter of nuclear matter, gp, gn are the 
partial number densities of protons and neutrons. Note that the account of the normal FL 
amplitudes in the case of the effective Skyrme interaction, being quadratic in momenta, is 
reduced to the renormalization of free nucleon masses and chemical potentials. Expressions 
for the quantities ^oo;^03) which enter in Eqs. with regard for the explicit form 

of the amplitudes Uq, U2 [jl9[ read 

c - c - 

'500 — 7^ A^oo, '^oa — /^03, 

2moo 2mo3 

where the effective nucleon mass moo and effective isovector mass mos are defined by the 
formulas: 

■ T^[3ti+t2(5 + 4x2)], (10) 



2moo 2moo 16 

.2 



h ag, , ^ , 

-^[ti(l + 2xi)-t2(l + 2a:2) 



2mo3 16 

being the bare mass of a nucleon. The renormalized chemical potentials /Uoo,/Uo3 should 
be determined from Eqs. (H), (|^) and in the leading approximation on the ratios T /sp, ^/sp 
have the form 

1 1 

/^oo = -j{^J'v^ f^n), /^03 = 2^l^P~ '^P'" ^ 2m ' ^^^"^ 

where kpp^^ = (3vr^^p,„)^/^, mp and m„ are the proton and neutron effective masses, defined 
as 

2 11 2 11 



moo mp rrin nios nip m„ 

The critical temperature of transition to TS superfluid phase is found from Eq. (1), deter- 
mining the energy gap, in the linear on A approximation. Considering, that the interaction 
amplitude Vi is not equal to zero only in a narrow layer near the Fermi-surface, |^oo| ^ ^ 
(we shall set 9 = 0.1/ioo) and entering new variables x = ^oo/a'-oo, 9o = 9/fioQ, Tfj_ = T//ioo, 
we present this equation in the form 
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1 = ! - (12) 
xftanh^ - + tanh^ -|, 



/ N , "^00 /^03 

^03 I^OQ 



mooPF [Z 

T^F = ^^^^3 , Pf= y2mooAioo 

The results of numerical integration of Eq. (|12D are shown in Fig. |l|. For small values of 
asymmetry a there exist such regions of large and low densities of nuclear matter, for which 
we have two critical temperatures. When a increases, these regions begin to approach and 
at some value a = ac {etc ~ 0.071) it takes place contiguity of the regions, so that we have 
always only two critical temperatures (for the regions, where solutions exist). For a > ac 
the phase curves are separated from the density axis and turn into the closed oval curves. 
Under further increase of a the dimension of the oval curves is reduced and at some a = am 
the oval curves contract to a point {a^ ~ 0.092). For the values a > am the triplet-singlet 
superfluidity fails. Note, that our results concerning two-valued behavior of the critical 
temperature qualitatively agree with the results of Ref . , where pairing of nucleons 
with the Paris NN potential was considered. 

As follows from Eq. (^, an equation determining the energy gap at T = has the form 

^M--i,Ev.i.,.)^ (13) 

X {sgn (E(q) + eo3(q)) + sgn (E(q) - ^03(q))} 

Passing in Eq. (|13D to integration on a layer, we arrive at the equation for determining the 
dimensionless gap y = AsoZ/ioo: 

i=f r°^4^{i+Bg„(V?T7-^.-v.)} (14) 

4 J -00 y/x^ + y^'^ ^ "^03 

(for a > it holds true mos > 0, /io3 < 0). The contribution to the integral gives the domain 
on X, for which the function, standing as an argument of the function sgn, is positive. In 
particular, such values of the gap, density and asymmetry parameter of nuclear matter are 
possible, that this function has no roots with respect to x. In this case the whole domain 
from —^0 to 6o gives the contribution to the integral in Eq. ([I^) and we arrive at the equation 
of the BCS type at T = with the solution y = 6'o/(sinh(l/5f)). 

Let us first find the solutions of Eq. (0) in the case when the chemical potentials /loo, /^os 
are given in the main approximation on A/ep- The results of numerical integration of 
Eq. (0) are presented in Fig. ^ In the case of symmetric nuclear matter {a = 0) we obtain 
the phase curve with one-valued behavior of the gap. For small values of asymmetry a there 
exist such regions of large and low densities of nuclear matter (excluding some vicinity of 
the point ^? = 0), for which we have two values of the energy gap, where one of these values 
is the solution of the BCS type and practically coincides with the corresponding value of the 
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gap for the case a = 0. For a less than some a^-, the energy gap demonstrates double-valued 
behavior in the intervals {Qrain, Q'min) and {g'^ax^ Qmax), while in the interval {g'^i^, dmai) it 
has one-valued behavior (see Table | for the values of various boundary points q). When 
a increases, the regions with double-valued behavior of the gap begin to approach and at 
a = Oc (the same etc as for the phase curves TJ^q)) it takes place contiguity of the regions with 
two solutions. For a> two branches of the phase curves are separated from the density 
axis and combine to one curve, beginning and ending in some points of the phase curve with 
a = 0. In this case the energy gap differs from zero only in the interval {gmin, Qmax), where 
it has double-valued behavior. When a increases further, the boundary points of the phase 
curves move towards and at some a = a'^ (not equal to am for the phase curves Tc{g)) the 
branches of solutions contract to a point. The value a'^, determines the maximum value of 
the asymmetry parameter, when TS superfluidity exists at T = (a^ ~ 0.179). 

Let us now find the solutions of Eq. (13) while accounting for the influence of the finite 
size of the gap on the chemical potentials ^oo,fio3- The results of integration of Eq. (p!4D 
in this case are presented in Fig. ^. Here the solutions of Eqs. (^4|), (H), and obtained 
for different a, correspond to the different parts of the dome-shaped curve, contained by 
the dashed lines of the same type for a given a. One can see, that taking into account 
the feedback of the finite size of the gap through the normal distribution functions /oo, /o3 
in Eqs. (^, (^) leads to the qualitative change: instead of two-valued behavior of the 
gap we have universal one-valued behavior. The first solution of the BCS type, obtained 
in uncoupled calculation, remains practically unchanged in self-consistent treatment of the 
gap equation (|1^ and with sufficiently high accuracy equals to its value at a = in the self- 
consistent determination. The second solution in uncoupled scheme, to which corresponds 
the smaller gap width, under simultaneous iterations of Eqs. (||), (^, (|14|) tends to the 
first solution of the BCS type. Taking into account the finiteness of the gap results in 
the reduction of the threshold asymmetry, at which superfluidity disappears, to the value 
a'^ ~ 0.029. Thus, in spite of the smallness of the ratio A/ep {A/ep < 0.12 for all densities 
q), the backward influence of pairing correlations is significant. This is explained by the 
fact, that if the quantity A in Eqs. (H), (^ differs from zero, then the absolute value of the 
chemical potential yUos increases a few times as against its value at A = 0. The increase of 
l/iosi is equivalent to the increase of the effective shift between neutron and proton Fermi 
surfaces, that leads to significant reduction of the threshold asymmetry. As at a 7^ the gap 
is finite everywhere, superfluidity arises and disappears under changing density by means of a 
first order phase transition. In principle, this phase transition can be observed in laboratory 
conditions under the study of intermediate-energy heavy ion reactions. If we assume that 
the final stage of the reaction can be described as an expansion of a compound nucleon 
system ||^, formed in a heavy ion collision, then under lowering density this disassembling 
phase can undergo a first order phase transition in density from the normal to superfluid 
state. 



The temperature dependence of the energy gap was studied in Refs. |11,13| with the 
use of the bare interaction in the gap equation in the form of the Paris and Argonne V14 
potentials, respectively. Our results agree qualitatively with theirs at T = 0. However, in 
our calculations with the effective density-dependent NN interaction we obtain the g cip clS cl 
function of density (not at fixed density) and this allows us to find new important features 
in behavior of the energy gap. 
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In summary, we studied TS superfluidity of asymmetric nuclear matter in the FL model 
with density-dependent Skyrme effective interaction (SkP force). In the FL approach the 
normal and anomalous FL amplitudes are taken into account on an equal footing and this 
allows us to consider consistently within the framework of a phenomenological theory the 
influence of medium effects on superfluid properties of nuclear matter. It is shown, that for 
some values of density and asymmetry parameter of nuclear matter the critical temperature 
of a second order phase transition demonstrates double-valued behavior, that agrees with 
the results of previous studies. In the case when the chemical potentials /iocyUos (half of 
a sum and half of a difference of the proton and neutron chemical potentials, respectively) 
are determined in the approximation of ideal Fermi-gas, the energy gap demonstrates for 
some values of density and asymmetry parameter the double-valued behavior. If we con- 
sider the feedback of pairing correlations through the dependence of the normal distribution 
functions of nucleons from the energy gap, then the energy gap drastically changes its be- 
havior from two- valued to universal one- valued character. In spite of relative smallness of 
the ratio ^/sf-, taking into account of the finite size of the gap in chemical potentials leads 
to the significant increase of absolute value of /ios and, hence, to considerable reduction of 
the threshold asymmetry, at which superfluidity at T = disappears. In self-consistent 
determination the energy gap at T = as a function of density has a finite width and 
normal-to-superfluid and superfluid-to-normal phase transitions should appear as a first 
order phase transitions in density. Among the other problems we note here the study of 
multi-gap superfluidity ||20| in asymmetric nuclear matter. 
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FIGURES 

FIG. 1. Critical temperature as a function of density. 
FIG. 2. Energy gap as a function of density in uncoupled calculations. 
FIG. 3. Energy gap as a function of density in a self-consistent scheme. 
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TABLES 



TABLE L The values of the boundary points (in fm '^), determining the intervals of double- 



and one-valued behavior of the energy gap in uncoupled calculations. 
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Qmin 


6min 


6 max 


Qmax 
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